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Abstract 

The nonperturbative mechanism of baryon-antibaryon production 
due to double quark pair (qq)(qq) generation inside a hadron is consid- 
ered and the amplitude is calculated as matrix element of the vertex 
operator between initial and final hadron wave functions. The vertex 
operator is expressed solely in terms of first principle input: current 
quark masses, string tension a and a s . In contrast to meson-meson 
production via single pair generation, in baryon case a new entity 
appears in the vertex: the vacuum correlation length A, which was 
computed before through string tension a. As an application elec- 
troproduction of A C A C was calculated and an enhancement near 4.61 
GeV was found in agreement with recent experimental data. 

1 Introduction 



The baryon-antibaryon (BB) final states in hadron reactions are rather typ- 
ical phenomena, e.g. in charmonium decays pp, AA etc. channels are sig- 
nificant pp. In e + e~ collisions the BB final states are carefully studied and 
display in many cases (A+A~,AA,pp) a nontrivial behavior near the corre- 
sponding thresholds, see [2] for a review and references. In B decays the 
produced pp pairs were observed with near-threshold enhancements [3j. In 
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Figure 1: Double quark-pair creation at points x,y in the field of heavy 
quarks cc. 



this paper we consider a rather general type of reactions, when a quarkonia 
state (QQ) decays into BB, where B{B) contain quark Q (antiquark Q). From 
dynamical point of view, the simplest case is the OZI allowed decay of heavy 
quarkonium into BB pair of heavy- flavor baryons, e.g. ip{nS) — » A+A~, 
which was experimentally observed first in jl] at one energy, and measured 
in [5] in the mass interval [4.5 -j- 5.4] GeV/c 2 . For this type of reaction the 
creation of two light quark pairs is necessary and one could expect some sup- 
pression in this channel. However, experimentally the suppression is quite 
mild, as was discovered in the reaction e + e~ — > A+A~ in jH |5]. 

A peak at the A+A^ mass around 4.63 GeV/c 2 was found in [5], and 
the nature of this enhancement is still obscure, however different explana- 
tions were suggested |7]. A discussion of possible mechanisms of similar 
phenomena in BB produced in B meson decays, was given in [§J. Below 
we develop a systematic theory of BB production in OZI allowed processes, 
which is actually a theory of double string breaking with BB emission, as 
shown in Fig.l. As it will be seen, this theory is a one-step development of 
the general approach of string breaking, given in [9]. To simplify matter we 
consider first the case of heavy quarkonium, decaying into heavy-flavor BB 
pair. 

2 The formalism 

The initial state of our problem is the heavy QQ state, where Q and Q are 
connected by a string. We are looking for a process, where two light qq pairs 
are created in the field of QQ., and hence the basic vertex is the 4q operator 
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in the static QQ confining field. As in the case of one pair qq vertex, it 
is sufficient to consider the light quark Lagrangian in the field of the static 
antiquark Q and static quark Q. 

This situation is shown in Fig.l, where a pair uu is created at the point 
(x, X4) and dd at (y, 2/4) with time growing from left to right. The string 
junction trajectory is shown in Fig. 1 by dotted lines and the string junction 
positions at each moment of time is defined as the Torricelli points in the 
triangles formed by space positions of (cud) and (cud). 

It is important, that points (x, x 4 ) and (y, 2/4) will be shown to be close 
to each other, and the string junction and anti-string junction are generated 
at one point in their vicinity, which considerably facilitates the picture of BB 
creation (while the latter is rather complicated in a two-step BB production). 

We start with the partition function of a light quark in the field of external 
current of heavy quarks QQ. 

Z = J DAD^jD^jexp [-(S + S, + S mt + S Q + Sq)] , (1) 
S, = \j d*x(F; u )\ (2) 



Si = -i J d 4 xip f (d + m)i) f , (3) 

S int = - j d A x^gA a t a ^. (4) 

Here / is flavor index, Sq and Sq refer to action of external quark cur- 
rents, of (possibly high mass) quark Q and antiquark Q. 

We exploit the background formalism [10J to split gluon field into confin- 
ing background and perturbative gluon field a M 

A^ = B^ + a M . (5) 

As in [5] we shall use the simplest contour gauge [TT] to express in 
terms of field strength^ 

r 11 ■ 

B »( x ) = / cHp(u)F itt (u)dv,i, a 4 = 1, = — , (6) 

JC(x) Xi 



1 Since the whole construction of S e f / for quark q in the field of antiquark Q is gauge 
invariant, the final result does not depend on gauge [12], and the use of contour gauge is 
a matter of convenience. 
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and the contour C(x) is going from the point x = (x, x 4 ) to the point (0, x 4 ) 
on the world-line of Q and then along this world-line to x 4 = — oo. Note, 
that our final result f lTT]) . f|T2|) will be cast in the gauge invariant form, which 
is the same for all contours, connecting points x, y to the world lines of Q 
(or Q). The independence of the resulting asymptotic expressions from the 
form of contours is shown in Appendix 3 of [13]. 
Averaging over fields B^, (F^), one can write 

Z = J J DipD^exp [-(S x + S eff )} , (7) 

where S e ff was computed in [T2j-[13j. Keeping only quadratic correlators 
and colorelectric fields for simplicity, one obtains (for one flavor) 

S eff = -- J d 4 xd 4 yip{x)^ 4 [ilj{x)ip{y)}^ 4 ^{y)J{x,y) (8) 

where [iptp] implies color singlet combination, and J(x, y) is expressed via 
vacuum correlator of colorelectric fields, 

J(x,y) = ^-(A 4 (x)A 4 (y)) = l* dm F dv t D{u-v). (9) 
iV c Jo Jo 

Here D(w) is the np correlator, responsible for confinement [15], 

jM^W) = (.6*8^ - S w 5, k )D(u -v) + 0{D X ) (10) 

and we have omitted the (vector) contribution of the correlator TJi, con- 
taining perturbative gluon exchange and nonperturbative (np) corrections to 
it. 

The properties of the kernel J(x, y) have been studied in [12] [13], here we 
only mention the general form 

J(x,y) = xy/(x,y)e^-^ 2 / 4A2 J D(0), (11) 
where we assumed the Gaussian form for simplicity D(x) = D(0)e Ja 7 , and 

/(x,y) = £ ds £ dte-t"-*F , x,y = ^, ^, (12) 
at small x,y, /(0,0) = 1, while asymptotically 



where 9 is the angle between x and y. Note also, that -D(O) and A are 
connected to string tension a 

a = 2ttA 2 D(0) = - J D{x)d 2 x. (14) 

We now turn to the effective action (JHD, where we write explicitly all flavor 
and color indices. In the latter case one should carefully restore the gauge 
invariant combinations, derived in [12], using parallel transporters $(w,t>) = 
Pexp(/J A^dz^) and we denote 

ipa(x)ipa(y) = $ a (x)® ab (x, Q, y)ip b (y) (15) 

with 

®ab{%, Q, y) = $ac(x, x 4 ; 0, x 4 )$ cd (0, £ 4 ; 0, y 4 )*db(0, 2/4, y, 2/4), (16) 
where is at the position of Q. Thus (|S]) can be rewritten as 
1 r 

s eff = ~2j d4xd4 y^l( x )^{( x )^l(y)^i(y)J(x,y)- (17) 

We take now into account, that A ~ 0.1 fm [14], [15] is much smaller, than 
all hadron scales, and one can integrate in ( 1TT|) over <i(x 4 — 7/4), using the 
form ( TTTT) . yielding 

5.// « -/d d 3 x C / 3 y(4 / (^)74^(x))(^(z/)74€(y))^(xy)/(x,y) 

(18) 

where we have used ( fl4l) and defined /(x, y) = so that /(x, x) = |^|, 

at large |x|. 

To proceed to the practical calculations with the realistic baryon wave 
functions, it is convenient to go over from bispinor to 2 x 2 formalism, as it was 
done in [16] for qq vertices, see Appendix 2 of [IB] . [Note, that the relativistic 
formalism for the hadron decay, developed in [0], [T7], and adapted for the 
baryon-antibaryon case in Appendix below, accounts for the full relativistic 
structure of hadrons, and is exemplified in the factor 2/123, which is the ratio 
of the vertex Z n factors for all hadrons. Below we follow a much simpler 
derivation in terms of 2 x 2 formalism, exploited in |16j .] 

We now take into account as in Appendix 2 of [16] , that each bispinor ip of 
light quark in (|18|) obeys the Dirac one-body equation (ap + j3(m + U))ip = 
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(e — V)ib, where U is the scalar confining interaction, U(x) = er|x — Xq|, 
and V corresponds to perturbative gluon exchanges; therefore one can write 

ib = I I , where 

V w J 

1 . , 1 

w = 77 77 { <T P) V 777 77 T^P 17 (19) 

m + f/-y + £: V 7 m+(U-V + e) K ' 

where angular brackets imply averaged value for a given quark in the given 
hadron, in our case this refers to the average energy and potentials of a light 
quark in the produced heavy-light baryon , e.g. A c . We also introduce for 
antiquarks bispinors ib c and spinors v c ,w c ,ib c = (v c ,w c ). Therefore 

$ = C~ V c = ^C-y = V c 7274 ; 7, = -0on 

and 

= -i( v c ,w c )/3 1° o 2 ) ( I; ) = ~^ C '* C ) [v) = - l (v c w + w c v) 



where notation is used, v c 02 = v c , w c <j2 = —w c = — v c cr p" m+ ^y +£ ^ • Hence 



(20) 

TBI) can be written as (we omit below superscript c in spinors v c 
S eff = [ d^x^O^ 



(21) 

where 

1 , 1 _ <x(p + fr) _ <tP 

A = ; -crp + cr p ■ = = . 

m+(U -V + e) F ^m+^-V + e) 

(22) 

We now form the S'-wave baryon wave function, which can be written as a 
product of a symmetric coordinate part and antisymmetric spin-flavor-color 
factor A e |H 

ijk vo 

(23) 



2 We neglect the nonsymmetric coordinate part of wave function, which contributes less 
than one percent to the nucleon mass, see [18, 19, 20 for more details. See also [20] for 
estimates of fl in (l22l). 
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where abc are color indices, a/?7 spinor indices and fgh flavor indices. 

One can separate the cm. motion and define the set of bound state wave 
functions in the cm. system {^(xi — x 3 ,x 2 — x 3 )}, 



*(r rd) (xi,x 2 ,x 3 ,,t 4 ) = 6 Vnfori), (24) 

where £, 77 are Jacobi coordinates, which can be defined in the relativistic 



iEti— iPR, 



case as [18], (w< = (Jpf + mf), u + = J2i^i 



rj = (x2-x x )J — 7~~~~~~~ r, £ = «/ 2 , 3 , r(wix 1 + w 2 x 2 -(a;i + a; 2 )x3) 

(25) 

and *&„(£, 77) is expanded in the fast converging hyperspherical series, where 
the leading term (> 90% in the wave function normalization, see [18], [19] 
for details) is a function of hyperradius only, \l/ n (£, ry) ps ^(/o), where 

p 2 = E-(x,.-R) 2 = | 2 + r/ 2 . (26) 
i=i w + 

In what follows we shall be primarily interested in the charmed baryons, 
A c , H c , Q c and their orbital (and radial) excitations. As a first example 
we consider A c and take for simplicity only one (leading) component of wave 
function with singlet diquark made of u, d. The explicit forms of Aq for 
p, A, X}>S are given in Appendix 1. For A(A C , A&) one can write in obvious 
notation 

A\ c = N Ac -\=e abc c aa {i){(ud) - (du)) jkjbc (27) 
ijk vo 

where (ud) jkbc = Up h (j)dp c {k)ep, Si = -e_± = 1. 

As shown in Appendix 2, the gauge invariant matrix element in the cm. 
system of decaying charmonium state *$? ni (r) can be written as 

G(niPi,n a P 2 ,n 8 Ps) = (2n) 4 8^(P 1 - P 2 - P 3 ) J n fn 2 n 3 (p) (28) 
where J** 2 

■C 2 n 3 (p) = / Z/i23e ipr rf 3 (x - u)rf 3 (u - v)rf 3 (x - y)(* ni A** n2 * n3 ). (29) 



Here r = c(u — v), c = UJQ ^ +UJd , and \I/ n . are coordinate space spinor 
wave functions,while M. is defined as 

M = *(xy)f(x,y)K x K y (30) 

At this point one needs to calculate the matrix element of the operator 
K x K y between BE wavef unctions, which we write as 

(p^p ) 

(A B \K x K y (vQaiV Q )\A B ) = rj BQ ^ 2 V (va^a) (31) 

Explicit calculation yields coefficients t/bq, given in Appendix 1 for A c , Y^,P, S. 

It is more convenient to go over to momentum space in J ni n 2 n 3 {p)i and 
using Appendix 2, Eq. (1A2.16|) . one has (we omit the superscript red in 
(IA2.16I) here and in what follows) 



where 



jBB ( \ f - d^Px cPpy d 3 q x d 3 q y , . 

J ™(P) = J ^3 (27r)3(27r) 3 (27r)3(27r)3 <^P - P, " Py)>< 

x ^n 2 (P x ,Py)^n 3 (P x + q x ,Py + q y ) (32) 

y ™= 2V2N A n y M( ^ )riQA (33) 

and Ai(q x , q y ) is the Fourier transform of .M(x, y), Eq. fl30|) . modulo K x K y , 
the latter were taken into account in the prefactor of Ai in fl33|) . Also 
and v a are spinors for A~ and A+ respectively, while Oi refers to the spin of 
1 (QQ) n state. From (IA2.18I) one can write 



d d r 1 r 1 z^L^iSmL 

M(q x ,q y ) = - — —CT7C4\ 2 / / dsdt{2nf 5^ (tq x + sq y )e W . 
dq x dq y Jo Jo 

(34) 

Insertion of fl34l) into fl33l and fl32|) yields finally 



J Z 2n M = yfJ l Q dsdt J ^fe0 e - A2Q2 ^( C p-p,-p,)^(p^p,) 



(2tt) 3 (2tt) 3 (2tt) : 
*3(p s + sQ,P tf -*Q), (35) 
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where we have differentiated by parts in (q x q y )M(q x , q ), obtaining 

* = 4 ^r^ (ilk) ^ v ^- (36) 

Note the factor 4 in (|36|) . which comes from the accounting for two di- 
agrams in Fig.l, and two diagrams with interchanging u- and d- vertices 
between points x and y. 



3 Baryonic width of heavy quarkonium and 
the BE yield in e + e~ collisions 

Using J ni n 2 n 3 (p) in (|35|) . one can find the decay probability of the n\ state 
of QQ into BB in the states n 2 ,n3 respectively, 



dw(l 23) = | J*«^(p)| W* W (?i -V 2 - V 3 )^yr^ (37) 

where bar over |J(p)| 2 implies averaging over initial and sum over final spin 
projections; in our simple case \v^aiV\\ 2 = 1. We now turn to a more direct 
experimental process of BB production, namely e + e~ —> BB, which was 
observed in jUE]. The corresponding amplitude can be written as [2T] 



A n2n3 ( P , E)=J2 Cn(E)T nn2n3 = ]T C n (E) J*n 2 n 3 (P 
n n,m \E — E + w(E) J 



mn 2 n 3 1 

(38) 



Here E and w are matrices in indices n, m, of the QQ system, (E) r 



E S 



I T7i\ f y—^ Jnn 2 n 3 (p) Jmn 2 n- i (p) / Qn \ 

W nm (E) = / — — 2^ p ^ j—s. , (39) 

where n, m refer to the complete set of charmonium bound states, and 
J n n 2 n 3 (p) is overlap integral of the n-th charmonium state and n2,n 3 states 
of heavy-light mesons. In terms of A n2n3 the total crossection is 

Vn 2 n 3 (E) = J \A n2n3 (p,E)\\-^5(E - E n2n3 (p)) (40) 
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The factor c n (E) in (1551) accounts for the production of {QQ) n pair in the 
given process, in case of e + e~ — > (QQ) n one has 

47rae Q V6 
c ™ = ^ ^«(°)' 

and with the definition (l38l) 



67T ■ 12e 2 

Ai? n2n3 (£) = ^ 2 g El^(o)T nW2ra3 (E)|^r n2n3 , (4i) 

where 

dT n2n3 = kj^5(E - E mns (p)). (42) 
As a result, keeping only one state n in (1411) one has for 

94p(£w;(o) K„»i 2 . . 

AJ^CB) g (43) 

where J^ 27l3 (p) according to (IA2.26P can be written as 



jBB ?5/2jV4 aJ^)^ 
fi2 (^ + ^) 3/2 (l + 2«)3/ 2 ' 



(44) 



where parameters (3q, Ro, Y, C refer to (QQ) n and £>£> wave functions and are 
defined numerically in Appendix 2. 

The polynomial 1Z n (p) is due to (QQ) n SHO wave function, and is ob- 
tained in the way described in Eq. (A. 33). It can be approximated as 

n n {p) -2.1 ^1 - 0.034^ - 0.05 (j^J j (45) 

In (HU) C and T are values of C and T, (IA2.27p . flA2.28p averaged over 
(s, t) integration region. 

A rough estimate of Ai?gg in (j4"3~I) . using ( )45i) . near A+A~ threshold with 
4>(4:S) state for (QQ) n is 



Ai? (4) - » e- eXp( ~ 2 ' 5p2 ^ (46) 
^£|£-£ 4 + «; 44 (£)| 2 ' 
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with 

£ = 0.9 • 10 4 e^ 4 2 (0) fe) 2 /x2 (47) 

Taking the Breit-Wigner form near the mass E4 of ip(AS), one can write 
the cross section of cr(e + e~ — > A+A") 

a 4 (e+e- A+A;) = £ 2.2 • lO' 4 ^ eXp( ~ 2 V) r2 (mi), (48) 

where all energies are in GeV and a 4 in mh. 

To calculate £ in ( 1471) . one can use (22], [23] average energies in (A. 25) with 
u n <C co c , hence a ~ 1,6 ~ and average values of s,t, s = t « 0.5. This 
yields (5 ^ 0.375, T £ 0.3. 

Here i? 2 = 4(p 2 + M| ), and all momenta and energies are in GeV. One 
expects, that A = 0(1 GeV -1 ), as follows from the exponential fall-off of 
D(x) in [14J, this value of A can be varied for the Gaussian form used above. 
The values of fL = — . IT } m — , n = u,d, where averaging is done over total 

n m n + {U-V)+e n i ' ' to to 

baryonic state, and m u , rrid ~ 0, can be found from the analysis of baryons 
in [20], where for the light quark in a single orbital with o = 0.15 GeV 2 one 
has e U: d ~ 380 MeV, while (U — V) can be roughly estimated as 0.6 -r- 0.8 
GeV, which yields Q n « fi d « (1 1.2) GeV. Now for the estimate of ^(0) 
one can use calculations in [22], checked vs experiment, which give values of 
-Rns(O) = V^7n/Vis(0) with account of mixing with (n — 1) 3 -Di states. Masses 
E n { 3 Si), given in Table below, are calculated in [22] (upper line) using 
flattening, and without flattening in [23] (lower line). 



Table 



n 


1 


2 


3 


4 


5 


6 


E n , GeV 


3095 


3.682 


4.096 


4.426 


4.672 


4.828 




(3.068) 


(3.663) 


(4.099) 


(4.464) 


(4.792) 


(5.087) 


R n (0), GeV 3/2 


0.905 


0.735 


0.511 


0.459 


0.360 


< 0.445 



Inserting R n (0) = 0.46 GeV 3 / 2 for n = 4, from the Table one obtains a 
typical value of £ ~ 0.035A 4 , and the maximum of from ( )48l) is of the 
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order of 1 nb for A ~ 1 GeV -1 . This magnitude is in accord with average 
experimental data in [5]. 

Note also, that exp(— 2.5(p) 2 ) divided by the Breit-Wigner factor in f )48|) 
is a strong cutoff factor, which decreases by a factor of 2 for AE = 0.2 GeV 
from the threshold. As a result, one obtains from ( I48p the resonance-type 
behavior of with maximum around E = 4.61 GeV, and decreasing twice 
at E=4.7 GeV, as shown in Fig.2. This form and the magnitude of the cross 
section correspond to experimental data in [6]. Explicit calculations with 
realistic baryon wave functions are now in progress [23]. 



~ l 

X! 

a 




Figure 2: The cross section <r(e + e - — > A+A - ) in nb, estimated according 
to Eq. (f48l) with A = 1 GeV -1 as a function of total energy E in GeV (solid 
line), experimental points are from [5], dashed line is the best normalization 
fit of Eq. (HI with a factor of 0.52365. 



4 Summary and discussion 

We have formulated above the fully nonperturbative mechanism for BE pro- 
duction via double quark pair generation. This mechanism is an extrapola- 
tion of the meson-meson production mechanism by string breaking, studied 
recently in [9] . Our main motivation is to construct a nonperturbative theory 
of strong decays from the first principles, in a similar way, as it was done in 
the theory of hadron spectra in one-channel case, where all hadron masses 
are computed from the first principle input: current quark masses, string 
tension a and a s , see [221 [22] for recent results and references. The simple 
string-breaking mechanism was indeed established with the only parameter 
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a in [9], and appeared to be close to the well-known phenomenological 3 P 
model, in this way giving a theoretical foundation for the latter. In the 
present case of BE production, an additional (fundamental) parameter ap- 
peared: vacuum correlation length A, which is connected to the gluelump 
mass A = l/Mf pi], and the latter is again expressed via string tension 

a: Mgi = 6.15 y/a = 2.6 GeV. In this way our first principle program is 
supported, however the 0(A 4 ) dependence of cross-sections makes the theory 
very sensitive to a possible pro cess- depending renormalization of A. 

It is clear, that the same mechanism should work for the pair creation of 
other baryons, containing c(c) quarks, e.g. S C E C , S C S C , the only difference 
will be in coefficient £ and the dominant intermediate resonance ip(nS). In 
the general case one should sum up over n, as shown in (138|) and a complicated 
interference picture may appear. 

In this case, when only ^(4S I ) state was kept, and this state is not far from 
the A+A~ threshold, the resulting bump in Fig. 2 is rather prominent. In Fig. 
2 the predicted theoretical enhancement is compared with experimental data 
from [5]. One can see a reasonable agreement. 

The pair-creation mechanism, given in this paper, can be applied also 
to the case of light (strange) quarks (QQ). In particular, for the reaction 
e + e~ — > AA, EE, studied experimentally in [23], one can use the same equa- 
tion ( HSj) . where the role of the intermediate state can play 0(2170) and 
higher 0-mesons. 

Since radius of high-excited 0's is much larger, than that of ip(4:S), the 
corresponding \ip n (0) | 2 and (3 in (|47|) are smaller, and one expects the cros 
sections a(e + e~ — > AA, EE), to be order of magnitude smaller than those for 
A C A C , E C E C . This is supported by experiments in [25] . 

As for the case of the cross section cr(e + e~ — > A b A b , E5E5), our Eq. f f48]) 
applies here without modifications, except for the replacement of ip(4:S) by 
T(65*); the main suppression factor comes from E 3 in the denominator of ( 148]) 

and from f^H = |, while 1^(0) | 2 acquires factor 8.6, since Bs(0) = 1.35 

GeV 3 / 2 [26]. As a result the cross section for A^Aft production is one order of 
magnitude smaller than that for A C A C production. 

One should stress, that the theory, developed here and in [9], can be 
applied to string-breaking processes, where the energy transfer AE from 
"external" quarks QQ to the pair-production vertex is not large, AE ■ A < 1. 
In the opposite case one should take AE into account in the string profile 
function J(x,y) in fTTT]) . which strongly changes result, these effects are now 
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under investigation. 

The theory, proposed above, is purely nonperturbative and therefore quite 
different from the mostly perturbative approach, developed before for BE 
production (see [8], [27] for discussion and references). In this respect two 
approaches complement each other and the final goal can be to formulate 
the unified theory, where all particle yields are expressed via first principle 
constants. 

The author is grateful to A.M.Badalian, I.M.Narodetski, and M.A.Trusov 
for discussions and comments. Useful advices, suggestions and help of G.V.Pakh- 
lova in preparing Fig.2 of the present paper are gratefully acknowledged. The 
financial support of the Grant RFBR No 09-02-00620a is acknowledged. 
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Appendix 1 

Baryon total wave function in terms of individual quark spinors 



As = N B £ ^ £ C B {fgh\aMH{a, i)v$(b, j)v h 7 (c, k). (Al.l) 

abc vD Uk 



&abc 
abc i j k 

We use notations vf(a, i) = u + (a,i), vf(b,j) = d + (b,j), etc. Here 

2 2 

z, j, k = 1, 2, 3 and X^fc denotes permutations of 1,2,3, we also require j < k, 
then the proton wave function with spin up can be written as (color indices 
are suppressed for simplicity). 



I a - \ " &abc 



^pE^E^WIWM*)) - «j)d(k))}- (Al-2) 

abc VD ij k 

Here notation is used: (du) = d-U + — d + u^; N p = 
For A( J = ~) hyperon with spin up one can write 

A A = N A J2^J2s + (i)[(<j)d(k)) ~ (dUHk))], (A1.3) 



abc » ijA; 



and A/a = 77=. 



For E hyperons with spin up 
A s0 = AT E0 DD^KWWo + Ho) - 2«_(<)H ++ }, (A1.4) 

where A^o = | and [itcf]o = u+(j)d-(k) + u~(j)d + (k), 
M++ = u +(j) d +(k) + d+{j)u+(k). 



A s+ = A/ s+ 2 E ^{*+(<)([«**]o " 2 fl _(i)u+(;)«+(*)}, (Al-5) 

where A/ s + = y/2/6. For E~ one replaces in f ]A1.5j) all u quarks by d quarks. 
For S° hyperon one has 

A s o = Ax+(u <-» s). 
As a result of calculations of t]qb one obtains 

1 4 1 / A \ 

VsA = 1, ^ss = - = r/ U p = -, r/ M A = ^=. (A1.6) 
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In the last two coefficients one must take into account the contribution 
of u quark in the vector meson (uu) n , which gets into the singlet pair (ud) 
or (du). This contribution is proportional to vgcrPv ■ Pi — > ^v^aiV^P 2 . In 
case of rjuA the isosinglet component of {uu) n gives an extra factor of -4?. 



Appendix 2 

Relativistic derivation of the hadron — > BB amplitude 

We start with the fully relativistic formalism and we follow here the 
derivation given in [T7]. The initial stage is the point-to-point amplitude, 
which is the Green's function G^xy for cc state emitted at point 1 and 
baryons absorbed at points 2 and 3, while intermediate points x,y are the 
same as in the main text, i.e. where two light quark pairs of flavors / and g 
respectively are emitted, see Fig.l. 

One can write this amplitude as 

J G 123xy d A xd*y = J d 4 xd 4 ytr{T l S Q {l,2)T 2 S f q (2,x)S^2,y)x 

x T x M{x,y)r y S f q (x,3)S 9 q {y,3)T 3 S g (3A)) = (0\j Q (l)j B (2)Mj B {3)\0) 

(A2.1) 

where S^q are light (q) and heavy (Q) quark propagators, and Tj are vertices 
for given hadrons, e.g. Ti = 7^ for 1 state of charmonia etc., while T x = 
T y = 74. Finally, 

^(x,y) = a(xy)/(x,y) (A2.2) 

and one should integrate ( 1A2.1I) over d 4 xd 4 y. However, the physical am- 
plitude of a hadron decay into two hadrons A(riiPi, n 2 P2, n 3 P 3 ) should be 
obtained from Gi 2 3 Xy in two steps: 1) first one should go from coordinate 
points 1,2,3 to definite momentum states Pi,P 2 ,P3, and 2) one should go 
from point-to-point amplitude to hadron-to hadron amplitude, which is ob- 
tained by amputating in the matrix element (IA2.ip the pieces (0|jj|rijPj), 
which are proportional to hadron decay constant. E.g. for a vector meson 

(0\jl\n,P = 0)=e k] J^ri n) (A2.3) 
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Proceeding as in Appendix 2 of [T7], one arrives at the expression 

A(n 1 P 1 ;n a P 8) n 8 p s ) = {2it)W\V x - V 2 - ? 3 )fijP) (A2.4) 

where 



Cn 3 (P) = ^-/l/l23rf 3 (x- M )rf 3 ( M -^)rf 3 (x- 2/ )vI/ ni ( U -v) ( T(xy) e i P r /(xy)x 

x * n2 (x-u,y-u)* n3 (x- v,y- v), (A2.5) 

and r = c(u — v), c = - + ^ c + , where Ui is average kinetic energy of quark 
i in the hadron. Here \I/ ni are coordinate parts of wave functions^, and |7 12 3 
is computed as a ratio of total trace and hadron Z{ factors, (see Appendix 2 
of [TF] for details) 

2/123 = ^=^, Z x = tr{T x k Q T x k Q ) (A2.6) 

V -^1^2^3 



Z fc (A; = 2, 3) = tr(T k f[ ^—^±T k ) (A2.7) 

8=1 ^ 

^123^ = tr(r i A Q r 2 (A g -7 4 A g )(A g -7 4 A g )r 3 AQ) (A2.8) 

= m 9 - ip^ + ^74 = m,- - - ^74 . , 

9 2u q 2a;, ' 1 ' ' 



and 



Here cu g = (yrn^ + p 2 ), where the average is for the given hadron n. 
Examples of 2/123 for meson—)- 2 meson decay are given in [161 E] • 
A much simpler derivation can be made in the so-called Dirac formalism, 
introduced in [16]. In this case the final expressions are given in the form of 
2x2 matrices and it is convenient in this case to write in (1A2.4j) the Dirac- 



reduced expression J£^ 2n3 instead of J^ 2 „ 3 , and the former is best written 



3 We assume here for simplicity, that a relativistic state can be described by only one 
scalar function, otherwise one has to sum over all terms with coefficients 2/123 > s P ec tfic for 
each term i 
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in the momentum space (first in the simpler case, when A4(x, y) in (1A2.2D is 
taken as an effective constant M). 

y ' red (2-K) 3 (27r)^ " l( ^ CP ~ Px ~~ p y^ n '^ Px > Py)^ri 3 {-Px, ~P y ) 

(A2.10) 

where 

y' red = tr{V^Y^K{ Vx )MK{p y )Y^}, (A2.ll) 

and K defined in ([22]). 

Note, that flA2.1ip has the same structure, as Eq. (B3) in [16], with (crq) 
replaced by KMK in (IA2.11I) . 

As follows from the Table VII in [IB] , the reduced vertex = ^crj 

for 3 Si state (1 ) of charmonium, while for , k = 2, 3, one must choose 
the appropriate baryon vertex, which is given in Appendix 1. To illustrate 
our procedure, we consider a simplified example (where color indices are 
suppressed, but the final result coincides with the exact one) of a baryon, 
consisting of singlet (ud) pair plus c quark, as 

( r £ed )a/?7 = ) a'/3' 7 ' = ' (A2.12) 

Now, introducing P in K(P), so that K = with P = p + "p, 
eK{P x )K{P y )e = e^K Pa {P x ) 

one obtains 

(eKKs) = (A2.13) 

and finally 

^=-Tmk {x+, " x)M (A214) 

and for Q x Q y can be assigned the values Q u , Qd (or vice versa) with Q Uy d = 

m u>d + (U — V) + E u ,d- 

However, Q Uj d in (1A2.13j) can be easily extracted from K in Eq. ([22]) of the 
main text, and one can see, that in the approximation, when the denominator 
in K is kept constant (independent of x or y) K x ~ cr(p x +p x ) = 0. Therefore 
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one must now take into account the coordinate dependence of J\4(x,y) in 
(lA2.2j) and we write 

M(x,y) = M{*- u,y - u) = / M^, q„) |^ e ^(x-u) + iq y (y-u) 

(A2.15) 



and (IA2.10D . (IA2.11D are replaced by 

X ^ri2 (Pxj Py ) ^ri3 ( Px 9a;) Py ^y) 

where y' red in (1A2.10I) is replaced by y red = y red M(q x , q y ), and y red in (1A2.14I1 
Here M.(q X) q y ) ls the Fourier transform of A4(x, y) ( 1A2.2I) . 



M(q x ,q y ) 



d 3 xd 3 y - ^= [ [ dsdtexp 

2\y/TT Jo Jo 



AX 2 



e -iqx-iq y y_ 
(A2.17) 



Performing the integrals, one obtains 

^ q » H ^S (2A)3 n' 

(A2.18) 

Insertion of (1A2.18j) into (1A2.16I) yields (after integrating out 5- function in 



dsdt{27c) 3 5 {3) {tq x +sq y )e 



(s + t) 3 



(IA2.18I) and differentiating by parts) 



where y red is now 



d 3 p x d 3 p y d 3 Q 
(2vr) 3 (2tt) 3 (2tt) 3 

^(p^PyWPar + sQSPy-tQ), 

_ rcd = 3 • 2 1 / 2 A 2 7T r(x + ^ X ) a'Q 2 



*l(<P-P*-P V ) X 

(A2.19) 



(A2.20) 



To estimate the integral in (1A2.19j) we use the oscillator wave functions 
for ^i, and oscillator form of hyperspherical wave function for ^3, 
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^(x) = P(x)exp ( ) , * 2 (p) = iVexp ( , (A2.21) 



where p 2 = ^ + r7 2 = ELi(x«-R) 2 ^, w f = + m?), w + = ELi<* 

In p-space one can use (see (25) and [18] for relations between standard 
and Jacobi coordinates) 



v|> 2 = N 2 exp(-R 2 (pl + p 2 )), N 2 2 = (8irR 2 Y (A2.22) 

and P§ = \/W c (Vx +Py)' Pr? = 73(Py ~Px)- 

Relation between average (p 2 ) and i? 2 , is (p 2 ) = 6-R 2 , = (r|). For ^ one 



can first take for simplicity ^i(p) = iVi exp . 

Now one can integrate in (1A2.19j) over d 3 p x d 3 p y , which yields 

J(p,Q) = J fjt (^ ^(^ ~ P^ ~ P«)*2(P«> Pj^3(P !E + gQ, Py ~ tQ) 



3 



A^ 2 r= = , exp(-H), (A2.23) 



where 

\2 j2 /„ , ,\2 



A) 

and 



d 2 (s + ty 2 2 



+ afl 2 Q 2 (s 2 +t 2 )-2fe# 2 Q 2 st--f - v 7 Q 2 fl 2 (a-fr) (A2.24) 



{UJ \ +Uc \ b=^, d, = J- + Rl( a +b); d 2 = -^+Rl(a+b)Q(s-t). 

(A2.25) 



Finally the integration over d 3 Q can be done in ( 1A2.19I) . yielding 



(47r) b Jo Jo 



7i 



1 3/2 



(A 2 + CRl)diRl(a - b) 



ex P [-Rl(cp) 2 T] 
(A2.26) 



22 



where 



a + b ( a + b \ R, 



T l + 2(5 2 R 2 {a + b) [l + 2(3 2 R 2 (a + b)J 4(A 2 + CR 2 ) ' (A2 ' 27) 
and 

C-a(s +t)- 2bst - — ^(a - b) - 2{1 + 2p 2 R 2 (a + b) y ( A2 - 2 8) 

3-2^X Mx +a lX ) 
y N c Q u Q d y ' 

In ( IA2.26I) iVi = (|pr) 3/4 , if the SHO for ip x is used. 

For (n 3 Si)cc state in the same oscillator basis one should use instead, as 
in EL6] 



* 1 (n 3 S 1 ) = 4=R s n HO (p , P )- 1 ) " (27T):V " 



2>^l)! e -^ Ll / 2 /V 
4^" \T(n + l) e VA 2 



^0/ 
(A2.30) 

normalized as / |^i| 2 (^p- = 1- 

To understand the structure of the obtained result (1A2.25[) . one can use 
the limit of large mass tuq, i.e. u c ^> u n , n — u,d, which yields a = 1, b = 0. 
Another useful limit is a) 2P 2 R^ 3> 1, which is achieved for large size baryons 
as compared to the radius of charmonium: note, that Rq = \{r 2 B ) ) while 
(3q(2S) = 0.46 GeV [16J. In this case one obtains 

c ^ (s2 + (2) - T ^-whcWY (A2 ' 31) 

In the opposite limit: b) 2(3qR^ <C 1 one has 

c _ 2 + (2 _(^_(^^ r ^__3_ (A232) 

In both cases one can use our result (IA2.24I) for the simple Gaussian form 
of to derive the final result for a more complicated function ( 1A2.30I) . by 
using 
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or directly, introducing in (1A2.23P the ^ given in (IA2.30|) . 

With the simple exponential function one has Ni = (j^j and for 
/2q/?o «C 1 one has an estimate 

Finally, one should take into account also, that y Tcd = y^(Q), while the 
total coefficient should be 



v&i = vStiQ) + v£(Q) + v£(Q) + v£(Q) = 4y. 



du 



(A2.35) 



and this is the value, which should be introduced in (IA2.25j) instead of 



,.red „,red , ,-red 

v > y ^ y total- 
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